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Recently, it has been proposed that exotic one-dimensional phases can be realized by gapping out 
the edge states of a fractional topological insulator. The low-energy edge degrees of freedom are 
described by a chain of coupled parafermions. We introduce a classification scheme for the phases 
that can occur in parafermionic chains. We find that the parafermions support both topological 
phases as well as symmetry broken phases in which the parafermions condense. In the presence 
of additional symmetries, the phases form a non-Abelian group. As a concrete example of the 
classification, we consider the effective edge model for a v — 1/3 fractional topological insulator 
for which we calculate the entanglement spectra numerically and show that all possible predicted 
phases can be realized. 



One-dimensional quantum systems can support topo- 
logically non-trivial phasesP™ These are gapped phases 
whose ground states do not break any symmetry of the 
Hamiltonian, yet they cannot be smoothly connected to a 
trivial (site-factorizable) state. Some of these phases also 
possess fractionalized edge states, which cannot disap- 
pear unless the bulk gap closes. In many cases, the topo- 
logical phase retains its non-trivial nature only as long as 
certain symmetries are kept; a prime example of this phe- 
nomenon is the Haldane phase of integer spin chainsP^HB 
Intriguingly, paired fermions can support a phase which is 
not protected by any symmetry (besides the intrinsic par- 
ity symmetry of fer mion s) , and is therefore robust to arbi- 
trary perturbations. 4 5 This is the one-dimensional topo- 
logical superconductor phase, which is characterized by 
fractionalized Majorana zero modes at its boundaries. 12 
These zero modes resemble the ones found in the 
cores of vortices in topological superconductors in two 
dimensions!-^ Other than the fundamental interest in 
this phase, it has been shown to be potentially useful for 
quantum information processing.^ It has recently been 
proposed that these phases can be realized in quantum 
wires proximity-coupled to superconductors EaEsu Signa- 
tures of Majoran a zero modes have been observed in re- 
cent experiments! 17 * 18 ^ 

These breakthroughs raise the question whether there 
is a richer variety of robust topological phases that can be 
realized in strongly-interacting one dimensional systems. 
Recently, it has been argued by Turner et alP, and simul- 
taneously by Fidkowski et alP that a one-dimensional 
system of fermions with arbitrary interactions and no 
symmetry (other than fermion parity symmetry) can sup- 
port only two phases: the trivial phase and the topolog- 
ical superconductor phase. 

In this paper, we discuss a possible way around this 
no-go argument. One starts from an effectively one- 
dimensional system which lives on the edge of a higher- 
dimensional fractionalized phase (e.g., a fractional quan- 
tum Hall phase). Therefore, the elementary degrees 
of freedom in the one-dimensional system need not be 



fermions or bosons; they can be fractionalized anyonic 
quasiparticles of the underlying two-dimensional phase. 
This leads to new classes of topological phases which can 
never be realized in a strictly one-dimensional system in 
the absence of symmetry. We focus on the recently dis- 
covered example of counter-propagating fractional quan- 
tum Hall edge states which are gapped by proximity 
coupling to a superconductorfSaEU This system can be 
mapped to a one-dimensional chain of particles with 
parafermionic statistics 2 ^!, of the type discussed recently 
by FendlejJ^U. This proposal generalizes an earlier setup 
that consists of edges of a quantum spin Hall insulator 
coupled to a superconductor! 2 ^ The boundary between 
regions of the edge in different topological phases support 
new types of zero modes with non-abelian properties.^ 

In this paper we classify the possible phases of this sys- 
tem on general grounds, and demonstrate the stability 
of the zero modes that occur between different phases. 
The different phases can be distinguished by the sym- 
metry properties of their entanglement states, which are 
reflected in characteristic degeneracies of their entangle- 
ment spectra. 

This paper is organized as follows: In Sec. [I] we re- 
view the arguments leading to the classification of gapped 
fermionic phases in one dimension. In Sec. [IT] we discuss 
a one-dimensional model of parafermions and introduce 
a classification scheme for parafermions in terms of the 
possible fractionalizations of the symmetry operator. We 
derive an addition rule for combining chains in differ- 
ent phases, and discuss how the classification is enriched 
when additional symmetries are present. In Sec. |III| we 
first review how to generate parafermions on the edge 
of a 2D fractional topological insulator. We then intro- 
duce a setup which allows to realize different topological 
phases on the edge of a v = | fractional topological in- 
sulator and study the effective edge model numerically. 
The different phases are identified by their characteristic 
entanglement spectra. We summarize the main results 
and conclude in Sec. IIVI 
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Figure 1: Illustration of a bipartition of a chain into a segment 
(S) of length Ls and an environment (E). A symmetry op- 
eration Q acting on the important eigenstates of the reduced 
density matrix fractionalizes into two operators Q A and 51 s . 



I. CLASSIFICATION OF SYMMETRY 
PROTECTED TOPOLOGICAL PHASES 

We begin by reviewing briefly the arguments of Ref. 
leading to the classification of gapped phases of interact- 
ing fermions in one dimension. The main idea is to study 
the properties of the entanglement in the ground state 
wave function \ipo)- We introduce a bipartition of the 
system into a segment S of Ls consecutive sites and an 
environment E (see Fig. [TJ . We assume that L$ is much 
larger than the correlation length £. We will now study 
the eigenvalues pf and the corresponding eigenstates \<j) a ) 
of the reduced density matrix p s = Tr^l^o) (^ | for the 
ground state \ipo). The eigenvalues pf are directly related 
to the so-called entanglement spectrum^ ef = — logpf. 
As the system is gapped and has short range correlations, 
only a finite number of states \<j> a ) have an eigenvalue 
Pa > e for a ny finite e > 0. This is directly related to 
the area lawp3 In other words, only a finite number of 
eigenstates contribute significantly to the state |V>o)- We 
will refer to these states as the important states. It was 
realized in Ref. [5| that any symmetry operation f2 act- 
ing on the important states \4> a ) can be split into two 
parts ft ~ £l A Q B which are only acting on the left and 
right site of the segment. Different topological phases are 
then characterized by the algebras of the fractionalized 
symmetry operators Q A and f2 B . 

Let us consider the example of a one-dimensional spin- 
less superconductor with fermion parity symmetry de- 
scribed by the following Hamiltonian 



the Hamiltonian Eq. ([I]) can be mapped (up to an addi- 
tive constant) to 



H = 



L-l 



.7 = 1 



(4) 



The operators jj are Majorana operators which obey the 
relations {7i,7j} = 25,,-, 7. = 7J. We can identify two 
topological phases which are characterized by Q either 
being fractionalized into two bosonic or two fermionic 
operators Q A ,Q B acting on the edges of the segment, 
i.e., we find two phases with Q A Q B = e l ^Q B Q A where 
p = 0, 7r. For any system with p = it, all eigenvalues 
p a come in degenerate pairs. To see this, note that Q 
and the reduced density matrix ps can be diagonalized 
simultaneously. If p s \4> a ) = pfl^a) and Q\<f> a ) = q\4>a), 
then the state \4> a ) = Q A \4>a) has the same eigenvalue 
for the reduced density matrix p \ip a ) — Pa\ipa) but 
QlV'a) = — Qxl'ipa) and thus \ip a ) is an orthogonal eigen- 
state. Consider now the Hamiltonian in Eq. ([!]) with 
u = 0, t = 1. In this case, it is rather easy to show 
that there are only two eigenstates \<f> a ) on the segment 
S (with a nonzero weight in the density matrix). The 
fermion parity of these states is given by Q = —iQ A Q B , 
where Q A = 71 and Q B = 72L and thus p = it as the two 
operators anti-commute. In the case u = 1, t = the 
ground state is a simple product state and we find p = 0. 
In the presence of both parity Q 2 — 1 and time-reversal 
symmetry T 2 = 1, eight different topological phases are 
found (see Ref. El for details) . 



II. PARAFERMIONS ON A CHAIN 

In this section, we generalize the concepts introduced 
in the previous sections to a chain of parafermions. Such 
a system was recently investigated by Fendley in Ref. [23] 
where it was shown that it exhibits exact zero modes, 
analoguous to the Majorana edge modes. 



Hn = 



' c ] c j+i 



H.c 



c h 



(i) 



with t,u > 0. The operator Cj (Cj) creates (annihilates) a 
spinless fermion on site j. The first term comprises hop- 
ping of fermions as well as the creation and annihilation 
of pairs of fermions while the second term acts as a chem- 
ical potential. The fermion parity operator Q = e i7r nj 



A. Parafermionic model 

We consider a one-dimensional model of Zjv 
parafermions with the following commutation relations 
for % < j: 



XiXi = e^XiXi 



(with 



) commutes with H because the total Furthermore, we have 



number of fermions iVtotai modulo two is conserved, 
defining the transformations 

72J-1 = Cj+c] 

l2j = -i(cj -ct), 



By 

(2) 
(3) 



Xf = 1. 



(5) 



(6) 



The parafermions are coupled between neighboring 
sites and we choose alternating coupling parameters t 



3 



and u. The Hamiltonian for a chain of 2L sites with 
open boundary conditions reads 2 ^ 



L-l 



3=1 



+ "E {^ /N ^3-1X13 + e-" /JV xWi) ■ (7) 
3=1 

This Hamiltonian conserves the charge modulo N, i.e.. 
it commutes with the operator 



Q = n (-e-^xh-^ 



(8) 



3=1 



which measures the total Zjv charge. 

It is useful to map Eq. ^ to the TV-state quantum 
clock model using a Jordan- Wigner transformation as 
shown in Ref. [231 



-e^lII^K ( 9 ) 

fe<j 




with the matrices 

/0 1 
1 




V 1 



1 

0/ 



r = diag(l,e 2 "/ JV , 



e 4 ™^,...^- 2 ^- 1 ^) .(10) 



Note that the operators Tj, Oj are defined on a lattice 
with L sites. The Tj, <7j operators satisfy 



of = rf = l 
of"* - 



The commutation relation between cr. r is given by 



(11) 

(12) 
(13) 

(14) 



It can be easily verified that the operators defined in 
Eq. (|9| fulfill the parafermionic algebra. Inserting the 
transformed operators in Eq. ^ yields directly 

L-l L 

H = - t J2 {'fa* + 'Ui'j) - « E { t 3 + t} ) , (15) 
3=1 3=1 

In the case of periodic boundary conditions, one has to 
pay attention to a phase change in the coupling of er^ 
and a\ which becomes 



X2LX1 



n 

,k<L 



<T\. 



(16) 
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Figure 2: Two different dimerization patterns corresponding 
to the two phases of the parafermionic chain. The trivial case 
(u = 1, t = 0) is shown in (a) where pairs of para fermions 
Xj are localized on each site, yielding a unique and gapped 
ground state. Panel (b) shows a non-trivial state (u = 0, t = 
1) with uncoupled para fermions at the two ends of the chain. 
This state is gapped in the bulk but has gapless excitations 
resulting from the unpaired para fermions and the edges. 



The phase diagram of the iV-state quantum clock 
model has already been investigated.^ It is in the same 
universality class as the Z n Villain model which was ex- 
amined by Elitzur et al. 29 . For N < 5, it exhibits two 
phases. For t > u, we obtain an ordered phase with 
an TV-fold degenerate ground state. If t < u, we are in 
the paramagnetic phase and the ground state is unique. 
The phases are separated by a critical point at t = u. For 
N > 5, however, this critical point is extended and a crit- 
ical phase emerges in between the two phases in a finite 
parameter region. The phase transitions into that critical 
phase are of Berezinskii-Kosterlitz-Thouless (BKT) type 
with an essential singularity in the correlation length. 
The phase itself is a BKT critical phase. 

To understand the phase diagram of Hamiltonian 
Eq. ([7]) , we follow Ref. [23] by considering two limiting 
cases. The case t — 0,u = 1 (Fig. [2|(a) has a unique 
and gapped ground state which is a factorized state 
where the parafermions \j form localized pairs on each 
site. The case t = l,u = (Fig. [2jb) is more interest- 
ing. Here, the parafermions form pairs between neigh- 
boring sites, leaving behind two unpaired parafermions. 
In the case of open boundary conditions, the unpaired 
parafermions reside at the ends of the chain and yield 
an TV-fold ground state degeneracy while the bulk re- 
mains gapped. This can be directly compared to the case 
of the fermionic model discussed in the previous section 
where unpaired fermionic modes appear at the boundary. 
In this limit, the individual terms of H commute with 
each other and can be diagonalized simultaneously. Any 

ground state |^o> satisfies -e nl/N x{jX 2 j+i Ho) = IV'o)- 
This implies that projection of the charge operator Q in 
Eq. Q onto the ground state manifold can be written 

as Q e g oc XiX2L- l n thi s sense, one can say that in the 
low-energy manifold, the charge operator "fractionalizes" 
into two a product of two operators localized at either 
end of the chain. Since Xi an d Xil both commute with 
the Hamiltonian, but do not commute with each other, 
the ground state is multiply degenerate. One can di- 
agonalize H and Q simultaneously, in which case there 
are N orthogonal ground states with distinct Q eigenval- 
ues. Acting with either x\ or X2L transforms the ground 
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Figure 3: (a) Overlap of a ground state of the Z4 
parafermionic chain (Hamiltonian |7| with N = 4 for dif- 
ferent L) transformed by \i an d the ground state of the next 
q-sector. The state |-0o) is the ground state for q = 1 and \tpo) 
for q — —i. (b) Matrix elements of xl between two degenerate 
eigenstates \<fi a ) and \4> a <) of the reduced density matrix p s 
which lie in g-sectors that differ exactly by one charge. Cal- 
culations are performed for different chains of length 2L with 
perodic boundary conditions. The segment S always includes 
the sites 1 to L. 



states into another, as can be seen from the fact that 
XiQ = e~ 2lT% t N Qx\ an d similarly for X2L- 

The decomposition Q e g oc X1X2L involving only oper- 
ators of site 1 and site 2L is only exact in the limiting 
case u = 0, t = 1. However, we expect Q c g to be ap- 
proximately given by a product of two operators which 
are local at either ends of the chain even when u =^ 0, 
as long as one remains in the same phase. Each of these 
edge operators is localized over a region which extends to 
a distance of the order of one correlation length from the 
boundary of the system. As long as the correlation length 
is finite, the exact edge operators at the left and the right 
ends of the chain have a finite overlap with xu X2L- To 
demonstrate this, we have diagonalized the Hamiltonian 



Eq. (15 I with open boundary conditions exactly for dif- 
ferent chain lengths and various values of {u, t}, and cal- 
culated the matrix element of x\ between ground states 
with different eigenvalues of Q . The results are shown 
in Fig. [3^ a). For t = 1, u = 0, the matrix element of x\ 



between ground states with neighboring eigenvalues of Q 
is unity. For u > 0, the matrix element decreases, but 
remains finite. The data clearly shows that the overlap 
does not depend on the system size but on the correlation 
length. As soon as the latter is comparable to the former 
(/3 s» 0.4 for L = 4), the matrix elements differ for larger 
system sizes. 

Instead of looking at physical edge states of a chain 
with open boundary conditions, we may also look at the 
"virtual" edges states we get by diagonalising the reduced 
density matrix ps- For this we consider again a bipar- 
tition of a periodic chain into a segment (S) of length 
Ls and an environment (E) as discussed in Sec. |lj The 
reduced density matrix p s commutes with the charge op- 
erator Q, and thus the eigenstates \<j> a ) of the segment 
have a defined charge q modulo N. As in the case of 
fermions in Sec. [I] we claim that the charge operator act- 
ing on the important eigenstates fractionalizes into two 
operators acting on the two ends of the segment. These 
localized operators then transform the important states 
with different charge into one another. To support this 
claim, we calculated numerically the matrix elements of 
x\ between important eigenstates of ps with g-eigenvalue 
that differ by a charge of one parafermion. The results 
are shown in Fig. [3j^b) . Again, the overlap does not de- 
pend on the size of the segment as long as the correlation 
length is smaller than Ls- This strongly supports the 
fractionalization of Q into a product of operators that 
are localized on either ends of the segment . 



B. Classifying phases by symmetry properties of 
the entanglement 

We now discuss the classification of the possible phases 
of a Zjv parafermion chain, generalizing the technique 
used to classify the phases of a fermionic chain (Sec. j|. 
These phases are best understood by first performing 
the Jordan- Wigner transformation Eq. (pi, mapping the 
problem to a Z^v spin chain. In the Ising/Majorana chain 
(corresponding to N = 2) , the symmetry broken phase of 
the spin chain corresponds to a topologically non-trivial 
phase in the fermionic variables, whereas the I12 sym- 
metric phase of the spin chain maps to the trivial phase. 

In the iV-state clock model, the Zjv symmetry can 
spontaneously break down to any subgroup of Z^r. The 
number of subgroups is the number of divisors of N. If 
N = p^ 1 . . .p[ l is the prime factorization of N, the num- 
ber of divisors is P = (Vj + 1). Therefore, we expect 
to find P distinct phases. It turns out that these phases 
fall into two distinct classes. One class is characterized 
by the presence of "fractionalized" edge modes (analo- 
gously to the non-trivial phase of the Majorana chain). 
These phases lack any local order parameter, and will 
be referred to as "topological phases". The other class 
can be described by the presence of a spontaneously bro- 
ken symmetry, even in the parafermionic variables. The 
two classes of phases can be distinguished by consider- 
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ing the ground state degeneracy on a chain with periodic 
boundary conditions: in a topological phase, the ground 
state is unique (as in the Majorana chain), whereas in a 
broken-symmetry phase it is multiply degenerate. 

In order to find all the possible topological phases, let 
us consider a chain of length 2L of ZAr-parafermions with 

a Hamiltonian H ^Xi> • • • > X2iJ that conserves the Zjv 
charge Q = Y[f=i ' —e~ ni / N xt;_i Yo „- I , hence 



and therefore 



[H,Q] = 0. 



(17) 



We now consider a bipartition of the chain with periodic 
boundary conditions into a segment (S) of length Lg with 
Ls <C 2L and an environment (E) as shown in Fig. [T] 

Just as in the case of fermions in Sec. [T] and as sup- 
ported by the numerical data in Fig. |3jb), we can write 
the operator Q acting on the important states \4> a ) as 
Q oc Q A Q B where Q A and Q B are local operators act- 
ing on a finite region near the left and right edges of S, 
respectively. The operator Q itself is a charge-neutral op- 
erator. However, Q A and Q B can carry non-zero charges, 
such that Q A Q — e 27rlk / N QQ A where k is an integer, and 
Q B Q = er 2lxlk l N QQ B . Different values of k correspond 
to different topological phases. Below, we derive a con- 
sistency condition which constrains the allowed values of 
k. 

We consider powers of the operators \ defined by 



X ni = X N/P '* 



(18) 



iei,...,l, 

and the corresponding charge operators 

q^sqW, (19) 

The original charge operator Q can be written as a prod- 
uct 



(Q ni )° 



(Q 



ni \a t 



(20) 



with integer a\ , • • • , ai . We can again fractionalize Q r 
like the total charge operator as 



gn» = (Q^) A (Q^f 



(21) 



We now study the ways in which Q ni can fractionalize. 
Similar to above, while Q ni is charge neutral, the opera- 
tors {Q ni ) A , (Q 7li ) B can carry charge Q ni . We may write 
{Q ni ) A as 



(Q ni ) A oc bosonic operators x (x™ ; ) 



(22) 



and correspondingly, (Q ,li ) oc [\ r ' 



Commuting (Q ni ) A and (Q ni ) B leads to 



(Q ,H ) A (Q n >) B =e 



J ""? fe . 2 {Q n *) B (Q n ^) A . (23) 

A 



There is also another way to determine ki. Since (Q ni ) 
creates riiki parafermions, commuting (Q ni ) with Q 
yields 



(Q n *) A Q n * = e-^ n - fci g™- (Q ni ) J 



(24) 



(Q"T (Q n T = e-^ n i k * (Q n ') B (Q n >Y 



(25) 



In order for ( 23 1 to be consistent with (|25[) , we obtain a 



constraint on fc,;: 

,2 (A 2 - ki) _ n.h (ki - 1) 



N 



Pi' 



integer. 



(26) 



The only two solutions for \ ki (ki — 1) are ki — or 
ki = 1 since fcj and ki — 1 cannot be divisible by the 
same prime. Hence, each Q ni can be fractionalized in 
two ways. 

The two possibilities of fractionalizing the operators 
Q ni together with Eq. ( 20 1 lead then to 2 l possible values 
of k given by 



ciiTiiki mod N, 



(27) 



where each of the ki can be chosen as or 1. Since there 
is no number that is a common divisor of all summands 



in (27), each individual choice of ki corresponds to a dif- 



ferent k. 

We now derive the structure of the entanglement spec- 
trum in the different phases. The ground state wave func- 
tion j^o) has a well defined charge modulo N and thus 
[p s ,Q] = [p s ,Q A ] = [p s ,Q B ] - 0. Hence, Q and p s 
can be diagonalized simultaneously. For k ^ 0, the eigen- 
values of p s come in degenerate multipletts which can be 
seen as follows. 

If we have a state \<fi a } with p s \<j) a ) — p& \4> a ) an d 



Q\4>a) = <la\4>a), then \ip a ) = Q J 



p S \^ a ) 

Q\1>a) 



Pa\lpa), 



has the eigenvalues 

(28) 
(29) 



From (29 1 we can see that \ip a ) and \(j) a ) are orthogonal. 

Every (Q Ui ) can only change the charge of the state 
by e 2 ^ LrLiki = e 27! " lki ^ Pi ' . Since different have no com- 
mon divisor, we can act D — Y[i=i [^i (pP ~ 1) + 1] times 
with Q A on the state \(f> a ) until we obtain a state with 
the charge q again. We identify this state as \4> a )- The 
eigenvalues of ps are therefore at least D-fold degenerate. 

We have identified 2 l different topological phases - cor- 
responding to phases in the clock model where the sym- 
metry is broken down to some subgroup of Zjy. The 
number of subgroups is the number of divisors of N which 
is Jl!:=i( r i + !)• This number is greater than 2 l if any 
Ti > 1. What happened to the remaining phases in the 
parafermionic model? Let us consider the example of a 
4-state clock model. Three phases can be realized, one 
with a ground state having the full Z4 symmetry, another 
Z2-symmetric one and a completely symmetry broken 
one. The first and the third correspond to the trivial 
and the non-trivial phase in the parafermionic chain, re- 
spectively. The second one is a symmetry broken phase, 
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even in the parafermionic variables. This can be seen as 
follows. The squared parafermionic variables commute 
with each other, 



[xlx-] =o, 



(30) 



thus they act like bosons. Hence, we can have symme- 
try breaking and long-range order in the parafermionic 
variables: 



(xlx}) o. 



(31) 



The argument that allows us to fractionalize the charge 
operator and identify the topological phase is no longer 
valid in the presence of symmetry breaking. Note that a 
long range order in parafermionic variables is not possible 
if they obey a non-trivial exchange statistics. This is 
because in a system with a finite correlation length we 
have in the bulk 



(XiXi+i) -> (Xi) (Xi+i) for I -> oo. 



(32) 



If the operators do not commute, this expression has to 
go to zero. Otherwise it cannot be consistent with 



(XiXi+l) = e 2m/N (x*+/X*) 



(33) 



To support the claim of a symmetry broken phase 
in the parafermionic chain numerically we checked the 
ground state degeneracies for open and periodic bound- 
ary conditions numerically on the example of a Z 4 chain. 
Recall that in the topological phase, the ground state 
degeneracy found for open boundary should vanish when 
choosing periodic boundary conditions whereas it should 
be the same for both cases in the symmetry broken phase. 
The Hamiltonian 



L-l 



Hsb — 



e (*y (x2 J+ i) +(xi J+ i) (x 2J ) 



L 

£ 

3=1 



-1X2.7 



A 2 jX2j-l 



(34) 



has for example a two-fold degenerate ground state for 
open and periodic boundary conditions. This indicates 
the symmetry broken nature of the phase. The Hamilto- 
nian 



L-l 



top 



3=1 



(35) 



has a four-fold degenerate ground state for open bound- 
ary conditions, whereas it is unique for periodic boundary 
conditions. Thus the degeneracy depends on the topol- 
ogy of the system. 



C. Combination of chains 

We now turn to investigate the topological phases that 
result when we combine parafermionic chains. As long as 



the inter-chain coupling is not too strong, the combined 
system is in a topological phase that is determined by the 
phases of the constituent chains. We derive an "addition 
rule" for chains in different phases. This is a generaliza- 
tion of the addition rules for phases of fermionic chainsP 

Let us consider two chains with parafermionic op- 
erators 02j-i!02j for the first one and ^ij-Xi^ii fo r 
the second. The setup is shown schematically in 
Fig. Qa). Our first task is to combine the two chains 
into a single effective chain with parafermion operators 
X4j-3,X4j-2,X4j-i,X4j, see Fig.gb). 

Suppose that the Hamiltonians of the two chains, 
and H^, are given . In order to construct the Hamiltonian 
of the combined chain, H x , we imagine the two chains 
embedded in a two-dimensional plane. The parafermions 
are thought of as abelian anyons which acquire a phase 
of e 27Tl / N when exchanged counterclockwise. We must 
then pay attention to the paths the anyons hop along. 
Hopping processes that go above or below the chain cor- 
respond, in general, to different operators. E.g., motion 
above the chain, and then motion back below the chain, 
corresponds to an anyon moving in a closed circle, and 
should therefore give a phase that depends on the total 
anyonic charge enclosed by the path. 



We define the operator 



At, 



to describe hopping be- 



tween sites i and j in the first chain, such that the anyons 
move below the chain. Similarly, ipjipj and x\Xj describe 
hopping processes between sites in the ip and \ chains 
that occur below the chains. After combining the two 
chains, however, we must consider the fact that the hop- 
ping process between cf> sites follows a trajectory that 
passes the tp sites from above (see Fig. Qc)). To take 
that into account, we define the following rule to map 
operators acting on the <j) chain to those of the x chain: 



(36) 



n e7ri/jv x 4fe -ixL i ■ 

i<k<j 



where have assumed that i < j, and that the hopping 
occurs between (f> sites with odd indices. In this equation, 

may replace {4i-i> X41-3} b y {02. > xL-ah respectively, 
corresponding to a hopping process originating from a <j) 
site with an even index, and similarly {foj—i, Xij-3} can 
be replaced by {</>2j, X4j-2}- 



The additional factor in brackets at the end of Eq. ( 36 1 



emerges as follows. The operator Xa%— 3X4^—3 
a process in which a particle hops below the combined 
chain. The fact that hopping between former cj) sites 
passes the ip sites from above is implemented by next 
hopping back below each i/j site which was passed, and 
then hopping forward again, this time above the ijj site. 
This operation leads to a phase factor of e - 47 ™/ w t the 
power of the number of parafermions at the tp site. This is 
exactly the factor in brackets in Eq. ( 36 1 . Note that this 
issue does not arise in Majorana chains (N = 2) where a 



describes 



7 



<hj-i 4>2 
O O 



o o 



(a) 



(b) 



o o 



lf>2j-l V2j 



o o „ * o o * x o o „ 



(c) o o 




o o 



In the last expression we omitted an overall phase factor 
that we get from commuting the charge operators. We 
obtain for the resulting state 



fctotai — fc</> ~\~ kjf, 2 k^k^p 



(40) 



The new fctotai fulfills the consistency condition Eq. ( 26 1 
since 



k 2 

"-total 



-fctotai = k^-k^ + k^-k^+A (k 2 - k^) {k. 



and each k^ 



kill) i 
(41) 

separately is a multiple of N. If 



N is a prime number, then there exist only two phases 
k = 0,1. Combining two chains with k<j, = k^ = 1, we 
find that the resulting phase is characterized by fctotai — 
0. Thus, the addition rule of phases has a Z2 group 
structure. 



Figure 4: Combination of two chains. The hopping in the 
two independent chains is shown in (a). After combining the 
chains (b) not all hopping processes occur below the chain. 
The anyons of former chain (j> pass above the sites of former 
chain ip (c). 



phase of either e +l7! or e~ l7r appears when hopping below 
or above the chain, respectively, so the additional phase 
factor always squares to unity. 

Following the same logic, hopping processes within the 
ijj chain do not need to be modified when mapping them 
to the combined chain: 



(37) 



A similar rule holds for ip sites with even indices. 

Let us now turn to study the addition rule of topolog- 
ical phases in the <fi and ip chains. The fractionalization 
rule for the charge operator of the cj> chain, Q(4>), can be 
affected by the considerations above. We need to map 
Q(4>) to the language of the combined chain \ operators 
according to the rule of Eq. ( 36 1 . If we name the \ °P~ 
erators on the former (j) sites as <// and the ones on the 
former if> sites ip', we obtain 



Q A {<!>)Q B {<t>) -+ Q A (4>')Q B {cj ) ')Q^') 



-2k d 



(38) 



where the original Q A ' B ((f>) operators are fractionalized 
according to Q A,B {<j)) ~ ±fc *, respectively. If we com- 
bine the two chains, we can calculate the k characterizing 
the fractionalization of the charge operator in the com- 
bined chain as follows: 

Qtotal = Q(tf>)QW 

= Q(4>')Q(i)')- 2k ^Q(ip') = Q^OW) 1-2 ^ 



oc Q\<i>')Q A W) l - 2k *Q B W)Q B W) 



f\l-2kd 



D. Presence of additional symmetries 

When there are extra symmetries in addition to the 
fundamental 7L?q symmetry related to charge conser- 
vation, a richer classification of phases arises. Inter- 
estingly, as we will demonstrate below, the phases of 
parafermion chains can have non-commutative group 
structures. This is in contrast to the symmetry-protected 
phases of fermions or bosons in one dimension which 
are always described by commutative groups, even if the 
physical symmetry group is non-commutative)^^ 

Consider, for example, Zjv parafermions, where N is 
prime, with an extra Z^ symmetry R that describes an 
internal degree of freedom (e.g., an orbital degrees of 
freedom). The phases of such a system are classified 
in terms of how Q and R fractionalize into terms that 
act on the edges of a segment: besides the parafermion 
charge of Qa, k, we also have to keep track of I, the 
parafermion charge of Ra- Altogether there are 2N pos- 
sible phases, since k is constrained to be or 1 (see 
Sec. II B), while I can be any integer modulo N. Note 



(39) 



that the second integer is similar to the one that classifies 
symmetry protect ed t opological phases of bosons with 
Zn x Zjv symmetry^! If there are two symmetries R, S, 
the charge of the fractionalized operators Ra under the 
action of Sa can be expressed as Sj^Rj^S^R^ — e 27Tll / N 
where I classifies the different phases. 

All the possible phases can be formed by combining 
copies of two systems, one in a phase that we label by X 
and the other in a phase labelled by Y. The X phase is 
of order 2, meaning that combining two systems in the 
X phase results in a system in the trivial phase. The 
Y phase is of order N. The phase X is the parafermion 
non-trivial phase we studied before, where all the parti- 
cles are assumed to have zero charge under R (i.e., the 
operators Ra/b into which R fractionalizes do not carry 
any charge); it is defined by k = 1 and 1 = 0. The other 
phase, Y, is defined by k — and I — 1, and it contains 
particles which have a charge under R as well as Q. We 



use the notation (k, I) to represent the combination of 
phases X k Y l . 

These phases follow the same multiplication rules as 
the dihedral group (the symmetries of an iV-sided poly- 
gon where X corresponds to a reflection and Y corre- 
sponds to a rotation). To show this, we must find how 
to classify a product of two states \<j>) ® \ip). Following 
the arguments of Sec. |II C| when the two chains are com- 
bined. R(4>) gets a string added onto it, given by a power 
of Q(^): R -> R A (cj)')R B {(t>')Q{ip'y 2l <t> (see Eq. @; 
recall that R A (<p) and R B {4>) carry a charge of ±f^, re- 
spectively). Thus, 

i?totai oc R A {<t ) ')R A {^)Q A {^)- 21 * x 
b')R B (ij/)Q B {ij/)- 21 * 



x R L 



(42) 
(43) 

(44) 



Together with the addition rules for k from Eq. ( 40 ) we 
find that 

(k$, 1$) x (fe^, lq) = (ktj, + k^ - 2k (t ,k t f,,l < p(l - 2k^) + l^). 

(45) 

In particular, we find for any prime number N that 
Y x X = X x Y^ 1 : the phases do not commute for 
N > 2. This shows that the group is the dihedral group 
(which is non-abelian for N > 2). 

As a concrete example, we obtain for the case of N = 3 
the following multiplication table: 



H4> 


(0,0) 


(0,1) 


(0,2) (1,0) 


(1,1) (1,2) 


(0,0) 


(0,0) 


(0,1) 


(0,2) 


(1,0) 


(1,1) 


(1,2) 


(0,1) 


(0,1) 


(0,2) 


(0,0) 


(1,1) 


(1,2) 


(1,0) 


(0,2) 


(0,2) 


(0,0) 


(0,1) 


(1,2) 


(1,0) 


(1,1) 


(1,0) 


(1,0) 


(1,2) 


(1,1) 


(0,0) 


(0,2) 


(0,1) 


(1,1) 


(1,1) 


(1,0) 


(1,2) 


(0,1) 


(0,0) 


(0,2) 


(1,2) 


(1,2) 


(1,1) 


(1,0) 


(0,2) 


(0,1) 


(0,0) 



This is exactly the multiplication table for the dihedral 
group of order 6. 



III. THE PHYSICAL REALIZATION OF 
DIFFERENT TOPOLOGICAL PHASES 

We begin by briefly reviewing a physical setup which 
gives rise to parafermionic quasiparticles at the edge of 
a fractional topological insulator. We then extend the 
setup by coupling of the localized quasiparticles through 
different media which allows to tune the system into dif- 
ferent topological phases. 




O Parafermions 
H FM 



SC 



Figure 5: Physical setup: The edge modes of the two- 
dimensional fractional topological insulator are gapped out 
by coupling to superconductiong and ferromagnetic domains. 
At the domain walls, localized parafermionic modes emerge. 



a stack of two FQH states, with filling fraction = 1/m 
for spin up electrons and — — 1/m for spin down elec- 
trons, respectively. On the edge, there are two counter- 
propagating edge modes. In analogy with the edge of a 
quantum spin Hall insulator^, these edge modes can be 
gapped by proximity to cither a superconductor or to a 
ferromagnet! 19 ! 21 !. The same physics can arise by bringing 
two counter-propagating FQH edge states close together 
and coupling t hem e ither directly or by proximity to a 
superconductor .1121221 p or concreteness, let us consider an 
FTI disk with an array of alternating L superconduct- 
ing and L ferromagnetic domains on its perimeter (see 
Fig.§. 

The low-energy theory of the edge in these setups has 
been shown to be closely related to the parafermionic 
chainPS At each interface between a ferromagnetic and 
a superconducting region, there is a localized zero mode. 
Zero modes of different interfaces satisfy "anyonic" com- 
mutation relations, equivalent to those of Eq. (|5| with 
N = 2m. It is therefore natural to identify each zero 
mode with a parafermionic operator. If two interfaces 
are sufficiently close to each other, fractional quasi- 
particles can tunnel between them. The resulting low- 
energy Hamiltonian is of the form of Eq. ([7]), where the 
"parafermion hopping" terms describe quasi-particle tun- 
neling processes between different interfaces. 

Below, we demonstrate for the example of N = 6 that 
an appropriately designed FTI system can realize all the 
phases of a Zjy parafermionic chain. 



Filling fraction v = 1/3 



A. Parafermions on the edge of a fractional 
topological insulator 

We consider the edge of a two-dimensional fractional 
topological insulator (FTIJpS which can be thought of as 



Let us consider the Zg-charge conserving model which 
arises at the boundary of an FTI with filling factor 
1/3. We now discuss the physical realization sketched 
in Fig. [6] which allows us to access four different topolog- 
ical phases. 
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The edge states of a v = 1/3 topological insulator are 
alternatingly coupled to superconducting and ferromag- 
netic regions as explained in the previous subsection. If 
the regions are short enough, fractional quasiparticles 
can tunnel between adjacent SC-FM interfaces. How- 
ever, since the interfaces are only separated by a small 
distance of vacuum, electrons can tunnel between them 
as well. In terms of the parafermionic degrees of free- 
dom this leads to a term coupling triples of parafcrmions 
on adjacent sites. By varying the length of the regions 
and the geometry of the edge, the amplitudes of the two 
types of tunneling processes can be tuned relative to each 
other. Putting the above tunneling terms together, we 
obtain the following effective Hamiltonian: 



a Z3 parafermionic chain and a Z2 (Majorana) chain. 
The corresponding conserved charges of these chains are 
given by the operators Q ni introduced in subsection |II B| 
The operator Q 3 is the conserved charge of the Majorana 
chain and Q 2 is the charge of the Z3 chain. 
Let us introduce the operators 



m = x 

7» 



Xi ) 



having the properties 



ViVj 

„2 



1, 



for i < j, 



Ho 



L-l 

E 



-iU 



(4) 

wi/e t 



3 / ^3 

X 2 j + 1 , 



xlj+i 



X 2 „ 



7i - 1. 
{7i»7?} = 2 % 
h,7j} = 0. 



(48) 
(49) 



(50) 
(51) 
(52) 
(53) 
(54) 



E h 



X23X23+1 + e m/6 xl J+ iX 2j 



pW»/6 t v 1 „-"/6 Y t v 
e A2j-1A27 ' e X27X27-I 



ZU 3 



3 



3 



(46) 



e 7ri/6 72 i 7 2j -+i^ i '72j+i 



h.c. 



(xtj-i) (x 2j ) -(xh) (x 2j -i) 



which describes the hopping of single parafermions or 
three parafermions at once. The operators xf obey the 
algebra of Majorana fermions: 



Hq may then be written as 

Ho = E [*i ( 

(i3 7 2j -72i+l+^3 72i-l72i) 



-2/: 



(55) 



To analyze this Hamiltonian, let us employ a "mean field" 
approximation in which we set 



{X?,X, 3 }=2<5, 



(47) 



72j72i+l r ?2j ? ?2j + l 



(56) 



From the classification described in the previous section, 
we know that there are four distinct phases in a chain 
of Iiq parafermions. Each phase is characterized by a 
specific degeneracy of the entanglement energy levels. We 
now show that all four phases can be realized by tuning 
the parameters in the Hamiltonian Hq (47 1. A chain 



P -W3„t 



of Iiq parafermions may be viewed as a combination of 



nmfin 



In this approximation we obtain an effective mean-field 
Hamiltonian Hmf = H 7 + H n with 



E 



i t 



^3 J 7 2j 7 2j+ i 



FTI v = 



+i ("3 - yA 2j _iJ 7 2 j-i7 2 j + h - c - . (57) 



FTI edge states gapped out by: 
h FM — SC 



and 



e tunneling 
quasi-particle tunneling 



H FM 

I I Vacuum 

Para- 
fermions 



^ = E 



^r 2i e-"/ 3 r,t 



+ yV 1 e-"/ 3 4 1 , 2j +h.c. 



(58) 



Figure 6: Physical realization of the Hamiltonian Ho ( 47 1 on 



the edge of a fractional topological insulator (FTI). Chang- 
ing the geometry can be used to tune the effective tunneling 
amplitudes for different particles. 



Hence, this effective model leads to two decoupled chains, 
a Majorana and a Z3 parafermion chain. Of course, this 
approximation is rather crude. However, we can use it 
to estimate in which parameter region the four phases 
occur. 
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(0,1,0,0) (0,1,3,0) (1,0,1,0) (1,0,0,0) 
(ti,u lt t 3 ,u 3 ) 

Figure 7: Entanglement spectra for different values of the 
parameters (ti, ui, <3, U3) in Hq. 



To prove the existence of the four phases for the Hamil- 



tonian Hq (47 1, we investigate the entanglement spectra 



for different parameter values numerically by exact diag- 
onalization of chains of length L = 6 in the clock model 
variables corresponding to 12 superconducting and ferro- 
magnetic domains, with open boundary conditions. The 
results are shown in Fig. [7J 

The expected degeneracies in the entanglement spec- 
trum are clearly reproduced. The ground states for 
(ti,ui,t 3 ,u 3 ) = (0,1,0,0) and (1,0,0,0) are fully dimer- 
ized states with the pairing of parafermions either within 
the same sites or between neighboring sites. The en- 
tanglement spectrum in these phases consists of a sin- 
gle multiplet of degeneracy 1 and 6, respectively. For 
(ii, Mi, £3, 1*3 ) = (0,1,3,0) we get that each level in 
the entanglement spectrum is two-fold degenerate, cor- 
responding to the non-trivial phase of the 7 (Majorana) 
sectoral, whereas the r\ sector is in the trivial phase. 



Similarly, for (tx, u\, ta, Us) = (0,1,0,1), each entangle- 
ment level is three-fold degenerate, corresponding to the 
non-trivial phase of the 77 sector. We have therefore 
shown that all of the phases predicted by our classifi- 
cation scheme are realized in the proposed model. 



IV. CONCLUSIONS 

In this paper we derived a systematic procedure for 
classifying the phases of a ID parafermionic system. The 
procedure is an extension of the previously introduced 
classification scheme of ID interacting fermions. We have 
shown that a ID parafermionic system supports differ- 
ent topological as well as symmetry broken phases. The 
topological phases are characterized by a fractionaliza- 
tion of the operator which measures the conserved charge 
on a segment. Different topological phases can be identi- 
fied by degeneracies in the entanglement spectrum. The 
symmetry broken phases are described by local order pa- 
rameters measuring the condensation of parafermions. 
We have derived an addition rule for chains in differ- 
ent phases. We furthermore showed that in the presence 
of additional symmetries the phases of parafermions can 
form non-commutative groups. As a concrete physical 
realization of the model we studied v — 1/3 fractional 
topological insulator which has been shown to yield a 
low-energy effective N = 6 parafermionic edge model. 
Here we showed that a simple physical setup with two 
both electron and quasi-particle tunneling allows us to 
realize all four possible topological phases. 
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